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f-->) Abstract. Most three dimensional constitutive relations that have been developed to describe the 

T— I behavior of bodies are correlated against one dimensional and two dimensional experiments. What is 

^^ usually lost sight of is the fact that infinity of such three dimensional models may be able to explain 
these experiments that are lower dimensional. Recently, the notion of maximization of the rate of 

^ entropy production has been used to obtain constitutive relations based on the choice of the stored 

^""5 energy and rate of entropy production, etc. In this paper we show different choices for the manner in 

^~~^ which the body stores energy and dissipates energy and satisfies the requirement of maximization of the 

^^ rate of entropy production that leads to many three dimensional models. All of these models, in one 

I— —I dimension, reduce to the model proposed by Burgers to describe the viscoelastic behavior of bodies. 
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Pr 1. Introduction 

t^ An observation of a phenomenon or a set of phenomena leads one to conjecture as to its cause 

Q and forms the basis for the crude first step in the development of a model. An experiment 
1-H is then deliberately and carefully fashioned to test and refine the conjecture. Unfortunately, 
;> this procedure is rendered most daunting as one is not usually accorded the luxury of being 
able to perform sufficiently general three dimensional experiments while the models that one 
would like to develop are fully three dimensional models. Most of the general three dimensional 
constitutive models that are being used in continuum mechanics have been developed on the basis 
of information gleaned from one or two dimensional special experiments. It does not take much 
mathematical acumen to recognize the dangers fraught in the process of such generalizations 
as infinity of three dimensional models could be capable of explaining the lower dimensional 
experimental data. Of course, one does not corroborate a three dimensional model by merely 
comparing against data from a single one dimensional experiment. One tests the model against 
several different experiments, but these experiments tend to be simple experiments in view of the 
extraordinary difficulties in developing an experimental program that can truly test the full three 
dimensionality of the model, especially when the response that is being described is complex. In 
order to obtain a meaningful three dimensional model on the basis of experimental data in lower 
dimensions, one needs to be guided by enormous physical insight and intuition. This is easier 
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said than done and in fact most models that are currently in use are based on flimsy and tenuous 
rationale. 

One might be tempted to think that the dictates of physics would greatly aid in the develop- 
ment of models from experimental data. For instance, the second law of thermodynamics could 
play a stringent role in the class of admissible models. Similarly, invariance requirements such 
as Galilean invariance could also provide restrictions on the class of allowable models. Unfortu- 
nately, the sieve provided by all such restrictions is far too coarse as it permits several models to 
go through that exhibit undesirable properties. 

While modeling, one might start directly by assuming a constitutive relation between the 
stress and other relevant quantities. This relation could be an explicit expression (function) for 
the stress in terms of kinematical variables as in the case of Hooke's law or the Navier-Stokes 
model, or it could be an implicit relation as in the case of many rate type non-Newtonian fluid 
models. Assuming such constitutive relations implies six scalar constitutive relations (in the case 
of the stress being symmetric). One could also assume forms for the manner in which energy is 
stored and entropy is produced by the body and determine the constitutive relation for the stress 
by appealing to a general thermodynamic framework that has been put in place (see the review 
articles by Rajagopal and Srinivasa [1], [2] for details of the framework). The framework casts the 
second law as an equation that defines the rate of entropy production (see Green and Naghdi [3] , 
Rajagopal and Srinivasa [4J) and appeals to the maximization of the rate of entropy production 
(while Ziegler [3] had appealed to such a requirement, the context within which he made such 
an appeal is different from that required by Rajagopal and Srinivasa [T], [2]). The general 
thermodynamic framework has been used to describe a plethora of disparate material response: 
viscoelasticity, inelasticity, twinning, phase transition in solids, behavior of single crystals super 
alloys, mixtures, inhomogeneous fluid, etc. While the method seems exceedingly powerful, there 
are some interesting nuances concerning its application that the modeler should be aware of, and 
in this paper by constructing explicit examples we illustrate these delicate issues. It is important 
to recognize that one can obtain the same constitutive relations for the stress by choosing different 
forms for the stored energy functions and the rate of entropy production (see Rao and Rajagopal 
[0] who develop the non-linear three dimensional Maxwell model by choosing two different sets 
of stored energy and rate of dissipation). In fact, it is possible that several sets of stored energy 
and rate of dissipation function can lead to the same model. We illustrate this by considering 
four different sets of stored energy and rate of dissipation to obtain the model developed by 
Burgers j7j , and these four choices are different from a previous choice made by Murali Krishnan 
and Rajagopal [8]. It is interesting to note that by making the choice of two scalar functions, we 
can arrive at a constitutive relation for the stress, a tensor with six scalar components. Many 
of the one- dimensional models that have been developed to describe the response of viscoelastic 
materials was by appealing to an analogy to mechanical systems of springs (means for storing 
energy), and dashpots (means for dissipating energy/ producing entropy), though in his seminal 
paper on viscoelasticity Maxwell ^ did not appeal to such an analogy. Within the context 
of these mechanical systems, it becomes clear how one can get the same form for the stress 
by choosing different stored energy and rate of entropy production functions as one can choose 
different networks of springs and dashpots to effect the same response. 

In 1934 Burgers [7j developed the following one-dimensional model by appealing to a mechanical 
analog: 

(1.1) a + pia + p2(T = qii + q2e, 



where pi and p2 are relaxation times, gi and g2 are viscosities, and a and e denote the stress and 
the hnearized strain respectively. A three dimensional generalization of that was provided by 
Murali Krishnan and Rajagopal [H], within the context of a thermodynamic basis that requires 
that among an admissible class of constitutive relations that which is selected is the one that 
maximizes the rate of entropy production. The second law merely requires that the entropy 
production be non-negative and one would expect the requirement of maximization of the rate of 
entropy to cull the class of rate of entropy production functions. As we shall restrict our analysis 
to a purely mechanical context, instead of making a choice for the rate of entropy production we 
shall make a choice for the rate of dissipation (the rate at which working is converted to heat) 
which is the only way in which entropy is produced within the context of interest. 

We shall assume that the class of bodies we are interested in modeling are viscoelastic fluids 
that are capable of instantaneous elastic response. A body that exists in a configuration Ht under 
the action of external stimuli, on the removal of the external stimuli could attain a configuration 
K,p(t), which is referred to as a natural configuration corresponding to the configuration Hf How- 
ever, more than one natural configuration could be associated with the configuration Kt based 
on how the external stimuli is removed, whether instantaneously, slowly, etc. The natural con- 
figuration that is accessed depends on the process class allowed. In this study, we shall assume 
the natural configuration that is achieved is that due to an instantaneous unloading to which the 
body responds in an instantaneous elastic manner. A detailed discussion of the role of natural 
configurations can be found in Rajagopal (TU] and the review article by Rajagopal and Srinivasa 
[T] . Even within the context of an instantaneous elastic unloading, it might be possible that the 
body could go to different natural configurations Kp.(t), i = I, ....,n. 

When one provides a spring-dashpot mechanical analogy for a viscoelastic material one obtains 
a constitutive equation that holds at a point, i.e., the point is capable of storing energies like 
the various springs and dissipate energy as the dashpots, but it also has to take into account 
the arrangement of the springs and the dashpots. The central idea of Mixture Theory is that 
the various constituents of the mixture co-exist. That is, in a homogenized sense at a point, the 
model has to reflect the combined storage of energies in the springs and the dissipation of the 
dashpots based on the way in which they are arranged. Of course, a point is a mathematical 
creation that does not exist, and what is being modeled is a sufficiently small chunk in the body. 
This chunk can store and dissipate energy in different ways. The point of importance is various 
arrangements of springs and dashpots can lead to the same net storage of energy of the springs 
and the dissipation by the dashpots. Put differently, the chunk can respond in an identical 
manner for different ways in which the springs and dashpots are put together. This is essentially 
the crux of the paper. We have five different three dimensional models, four that are developed in 
this paper and one that was developed by Murali Krishnan and Rajagopal [8] and all five three 
dimensional models could claim equal status as generalizations of the one dimensional model 
developed by Burgers. Recently, Malek and Rajagopal [TT] used the thermodynamic framework 
that we have discussed to obtain a model for two viscous liquids. In this paper, we have a more 
complicated mixture in that we have two different elastic solids coexisting with two different 
dissipative fluids. We do not allow for relative motion between the constituents, we assume they 
coexist and move together. 

The organization of the paper is as follows. In section 2, we introduce the kinematics that is 
necessary to the study and the basic balance laws for mass, linear and angular momentum. We 
also introduce the second law of thermodynamics. This introduction is followed by a discussion 



of four different models which all reduce to Burgers' one dimensional model in sections 3-6. We 
make a few final remarks in the last section. 

2. Preliminaries 

Let kr{B) and Kt{B) denote, respectively the reference configuration of the body, and the 
configuration of the body B at time t. Let X denote a typical point belonging to hr{B) and 
X the same material point at time t, belonging to Kt{B). Let Xkr denote a sufficiently smooth 
mapping that assigns to each X e hr{B), a point x G Ht{B), i.e., 

(2.1) x:=x«^(X,t). 

The velocity v, the velocity gradient L and the deformation gradient F^^ are defined through 

It immediately follows that 

(2.3) L = F«,F;^. 

We denote the symmetric part of the velocity gradient by D, i.e., 

(2.4) D := 1 (L + L^) . 

The left and right Cauchy-Green stretch tensors B^^ and C^^ are defined through 

(2-5) B«,^ := F«,^F^^, C^^ := F^^F«,^. 

Let K,p(^t) denote the preferred natural configuration associated with the configuration k^. We 
define F,^ as the mapping from the tangent space at a material point in Kp(t) to the tangent 
space at the same material point at Kt (see Fig. [I]). We then define 

(2-6) Bkp(j) := F«;^(^jF^^^^^, C^pj,, := F^^^^^F^pj^,. 

The mapping G is defined through (see Fig. [I]) 

(2- ' ) G := FK^^Kp(t) := F^^^^^F^^. 

We can then define the tensor Ck^-s^k^^jj in a manner analogous to C^^ through 

(2.8) C«;^_^«;pjjj := G G, 

and it follows that 

We shall also record balance of mass (assuming incompressibity) , balance of linear and angular 
momentum (in the absence of body couples): 

(2.10) divi^v) =0, pv = div{T^) + ph, T = T^, 

where p is the density, v is the velocity, T is the Cauchy stress tensor, b is the specific body 
force, div{.) is the divergence operator with respect to the current configuration and (.)^ denotes 
transpose. In addition, the local form of balance of energy is 

(2.11) pe = T.L~div{q)+pr, 




Figure 1 . Schematic of the natural configuration Kp(^t) corresponding to the cur- 
rent configuration Kt and the relevant mappings from the tangent spaces of the 
same material point in hr^ Kt and K.p(t)- 

where e denotes the specific internal energy, q denotes the heat fiux vector and r denotes the 
specific radiant heating. We shall invoke the second law of thermodynamics in the form of the 
reduced energy dissipation equation, for isothermal processes: 

(2.12) T.-D-pij\i,othermar.= ^>0, 

where ip is the specific Helmholtz potential, ^ denotes the rate of dissipation (specifically rate of 
entropy production). 
When one works with implicit constitutive models of the form 

(2.13) f(T,D)=0, 

where T is the Cauchy stress, or more general models of the form 



(2.14) 



V 

f I T,T, 



(n) 
V V 

T,D,D, 



0, 



(n) 

where the superscript V stands for the n Oldroyd derivatives [12j, and where T and D seem to 
have the same primacy in that the maximization could be with respect to T or D. However, the 
superficial assumption that T and D have the same primacy is incorrect as T (or the applied 
traction which leads to the stresses) causes the deformation (the appropriate kinematic tensor). 
In order to get sensible constitutive equations one ought to keep D fixed and vary T to find 
how a body responds to the stress that is a consequence of the applied traction. This comes 
up naturally in the development of implicit constitutive theories (see Rajagopal and Srinivasa 
[13], Rajagopal and Srinivasa [E]). More recently, Rajagopal [15] has discussed at length the 
implicit nature of constitutive relations. When one thinks explicity along classical terms of the 
stress being given explicity in terms of the kinematical variables, it is natural to hold T fixed 
and maximize with respect to the kinematical variable, in our case D. This is what is followed 
in this work. 



3. Model 1 

3.1. Preliminaries. Let k^ denote the undeformed reference configuration of the body. We 
shall assume that the body has associated with it two natural configurations i.e., configurations 
to which it can be instantaneously elastically unloaded and corresponds to two mechanisms for 
storing energy (within one dimensional mechanical analog - two springs). Interestingly, one can 
get from the reference configuration to the two evolving natural configurations denoted by Kp-(j), 
i = 1, 2 (see Fig. |2|, via two dissipative responses. Let Fj, i = 1, 2, 3, denote the gradients of the 
motioiiF from kr to Kp^(t), Hp^{t) to i^p2{t)-, and ^^^(t) to Kt respectively. Also, we shall define the 
left Cauchy-Green stretch tensors, 

(3.1) Bi:=F,Ff, ^ = 1,2,3, 

and the velocity gradients with their corresponding symmetric parts, 

(3.2) L, :=F,Fri, D, := ^ (L^ + Lf) , 2 = 1,2,3. 

Also, we note that [j 

(3.3) F = F3F2F1. 

Let us denote the gradient of the motion from Kp^{t) to Kt by F^; then, 

(3.4) Fp = F3F2, 
and 

(3.5) F = FpFi. 

The left Cauchy-Green stretch tensor, the velocity gradient with its symmetric part, correspond- 
ing to Fp are 

(3.6) B,:=FpFj, Lp := FpF;\ Dp := ^ (Lp + LJ) , 
respectively. 



Now, taking the time derivative of eqn. (3.5) we get: 

F = FpFi + FpFi 
(3.7) ^ LF = LpFpFi + FpLiFi 

^ L = Lp + FpLiF-i. 



Similarly, 


taking 


the time derivati 


ve of eqn. (3.4), we arrive 


at 


(3.8) 






Lp 


— L3 + F3L2F3 . 




Now, 






Bp = 


p p ^ p p 





(3.9) = FpFjLj + LpFpFj 

= BpLp + LpBp. 

In general, these are appropriate mappings of tangent spaces containing the same material point in different 
configm'ations. 
Henceforth, we shall denote F^j^ by F. 




>i(i) 



'P2(t) 



Figure 2. Schematic to illustrate the natural configurations for model 1. kr is 
the reference configuration, Hf denotes the current configuration, and i^pi{t), i^p2it) 
denote the two evolving natural configurations. The body dissipates energy like a 
viscous fluid as it moves from, kr to ^^^(j), and i^p^{t) to Hp2{t)- Also, as shown, the 
body stores energy during its motion from, Hp^it) to Ht, and Kpi(t) to Kj. 



Post-multiplying eqn. (3.7) by Bp, pre-multiplying the transpose of eqn. (3.7) by Bp, and adding, 
we obtain 

(3.10) 



V 

B. 



-2FpDiFj^ 



wher e Bp := Bp — BpL-^ — LBp is the Oldroyd derivative of Bp. In a similar fashion, using 



eqn. (3.8) and the relation B3 = B3L3 + L3B3, we get 



(3.11) 



B3= -2F3D2F[, 



where B3:= B3 — BsL^ — LpB3. This is same as the Oldroyd derivative of B3, when the natural 
configuration Kp^{t) is made the reference configuration. 
We also note from eqn. (3.3) that 



(3.12) 



and hence 



F — F3F2F1 + F3F2F1 + F3F2F1 



L = L3 + F3L2F3 ^ I -P t:^ T T7-1-I7-1 



F3F2LiF2-^F3 



I.B3 = I. (L3B3 + B3L[) 

= I. (LB3 — F3L2F3 — F3F2L1F2" F3 + B3L 



FsL^Fj 



F3F2-^LfF^F[) 



(3.13) 

= 2B3.D - 2C3.D2 - C3. (F2LiF2-^ + F2 ^Lf F^) . 

The relations derived, in this sub-section, are sufficient for the purpose of analysing this model. 
In the following sub-section, we shall constitutively specify the forms for storage and rate of 
dissipation functions, and then we shall maximize the rate of dissipation subject to appropriate 



constraints (incompressibility and the energy dissipation equation), to determine the constitutive 
relation. 

3.2. Constitutive assumptions. Let us assume the specific stored energy ip and the rate of 
dissipation ^ of the forr^ 



(3.14) 



ij = tP{B3,Bp), e = e(Di,D2) 



In particular, assuming that the instantaneous elastic responses from Hp-^^^t) and K,pj^{t) are 
isotropic, and in virtue of incompressibility of the body, we choose 



^iB,,B, 



/Us 



2p 



IB, 



3) + |(I.B, 



3), 



(3.15) 

and 

(3.16) e(Di, D2) = ViDi.Di + V2^2.'D2. 

The above assumption means that the body possesses instantaneous elastic response from the 



two evolving natural configurations (/«pi(t), /«p2(t)) to the current configuration Kt (Fig. pi); the 
body stores energy like a neo-Hookean solid during its motion, from Hp^^t) to Kt, and from Hp^^t) 
to Kf. In addition, the response is linear viscous fiuid-like, as the body moves from nji to Hpj^{t), 
and from one natural configuration (^^^(t)) to the other {Kp2{t))- 

Also, since we have assumed that the material's instantaneous elastic response is isotropic, we 
shall choose the configurations Kp^(t), Kp^(t) such that 



V, 



V. 



(3.17) i a — » 3, 3. p — , p, 

where V3, Vp are the right stretch tensors in the polar decomposition i.e., the natural configu- 
rations are appropriately rotated. 
Finally, using eqns. (3.13) and (3.17), we get 

(3.18) 

and similarly 



IB, 



2B, 



D - D2 - - (F2L1F2 1 + F^-^Lf F^) 



(3.19) 



IB. 



2Bp.(D 



Di) 



Substituting eqns. (3.15), (3.16) into (2.12) and using the relations in eqns. (3.18), (3.19), 
(3.20) 



T.D-/i3B3. 



1 



D - D2 - - (F2LiF2-^ + F^'^Lf F^) 



/ipBp.(D - Di) = T^iDi.Di + r/^D2.D2 



which on further simplification leads to 

(T - /iaBg - ^lpBp) .D + /i3B3.D2 + /XpBp.Di + ^B^. {Y^Li^^^ + F2-^L[F^) 

= r/;Di.Di+r7^D2.D2. 



One can also choose the rate of dissipation function to depend on the stretch i.e., of the form ^ = 
^(Di, D2, B3, Bp). The resuhing constitutive relations will be a variant of the relations obtained when ^ is 
of the form given in eqn. (3.141. The constitutive relations obtained by using ^(Di, D2, B3,Bp) have relax- 
ation times which depend on the stretch. Upon linearization, the two constitutive relations take the same form. 
Rajagopal and Srinivasa have discussed this issue for the Maxwell fluid in [16]. 



We shall assume that the body can undergo only isochoric motions and so 

(3.22) tr(D) = 0. 

Also, since the body can actually attain the two natural configurations, the incompressibility 
constraint implies that 

(3.23) tr(Di) = 0, tr(D2) = 0, 
where tr {.) is the trace of second order tensor. 



Since the right hand side of eqn. (3.21) does not depend on D, along with eqn. (3.22), we have 
(3.24) 



T=-pI + /isBg + flpB 



pi 



where —pi is the reaction stress due to the constraint of incompressibility. Hence, eqn. (3.21) 
reduces to 

(3.25) /isBs.Ds + /ipBp.Di + ^Bg. [Y^LiF^^ + F2 ^L[F^) = r/iDi.Di + r/^Ds.Da. 



Following Rajagopal and Srinivasa [16J, we maximize the rate of dissipation in eqn. (3.16) along 
with the constraints in eqn. (3.23), (3.25), by varying Di,D2 for fixed B2,B3. We maximize the 



auxiliary function $ defined by 
(3.26) 

$:=77;Di.Di + 77^2.02 



+ A2I.D1 + A3I.D2 



/^3i 



ViDi.Di + r7^D2.D2 - /X3B3.D2 - /XpBp.Di - ^Ba. {Y^U^^^ + Fa'^Lf F^) 



Now, setting 9$/9D2 = 0, t?$/t?Di = 0, and dividing the resulting equations by Ai and A2 
respectively, for Ai, A2 7^ 0, we get (also see appendix) 



A^sBg 



(3.27) 



-^— J 2,,D2 + -I, 



^^,B, + I (F^BsF,-^ + F,-^B3F2) = (^) 2v[B, + ^I. 

/i3B3.F2WiF2-' 



Using eqn. (3.27) in eqn. (|3.25|), we get 
(3.28) 



Ai + 1 1 



Ai 2 2r7;Di.Di + 2r/^D2.D2' 

where Wi := | (Li — Lf). Hence, 

T = -pi + /i3B3 + UpBp, 

(3.29) y (F^B3F2-^ + F2 IB3F2) + /ipBp = -p'l + rj.B,, 

fXsBs = -p"\ + ?72D2, 
where p', p" are the Lagrange multipliers with 



(3.30) 



/I //I 

-P = 3 [f^3tr{B3) + fiptr{Bp)] , -/ = -fi^triB^), 

^1 = 2(^)V„ ,2 = 2(^),^. 



Now, eqns. (3.10), (3.11) can be re- written as 

1. 



(3.31) 



Di = --v;1b,v;\ d, 



-y^' B3 V.:; 



Using eqn. (3.31 )& in eqn. (3.29 )c, and post-multiplying and pre-multiplying with V3, we have 

1 



(3.32) 



/isB^ 



-/i3tr(B3)B3 



V2 -X^ 
2 ^ 



and using eqn. (|3.4|), we get 
(3.33) 



In addition, using eqn. (3.31 )a in eqn. (3.29 )b, post-multiplying and pre-multiplying with V 



P' 



/U3 



BpBa + BsBp) + /ipB^ 



1 v^ ^ 

- [fistriB^] + /iptr(Bp)] B^ - ^ Bp 



Notice, from eqns. (3.32), (3.33), that the evolution of the natural configurations Hp-^a) and 



Kp2{t) are coupled. These two equations are to be solved simultaneously to determine their evo- 
lution. We shall denote fi^'B^, fipBp by Si, S2 respectively. Then, the final constitutive relations 
- eqns. (3.29 )a, (3.32), (3.33) - reduce to 

-pi + Si + S2, 



(3.34) 



T 

2 



1 

2 



S2S1 + S1S2) + S2 



1 7l2 ^P 

-tr(Si)Si- — Si, 
^[tr(Si)+tr(S2)]S2-|s2 



In the next sub-section, we shall show that the above constitutive model reduces to Burgers' 
model in one dimension. 

3.3. Reduction of the model to one dimensional Burgers' model. In this sub-section. 



we shall first linearize the constitutive model given by eqn. (3.29) (we shall use eqn. (3.29), here 



instead of eqn. (3.34), for the sake of simplicity) by assuming the elastic response is small (we 



shall define what we mean by small, precisely, later). Then we shall show that, in one dimension, 
the equations reduce to the one dimensional linear model due to Burgers (see eqn. (1.1)). 
Now, eqn. (3.29 )c can be re-written as 

^^^(Bs) - 



(3.35) 



/is (B3 - I) 



At3 



I + r/2D2. 



If the displacement gradient with elastic response is small, i.e. 

au(X,t) 



(3.36) 

then 

(3.37) 
and hence 
(3.38) 



max 
xgb 
ten 



IB, -II 



dX 



0(7), 7<1, 



0(7), 7<1, 



3,P, 



tr(B,) = 3 + 0(7' 



3,P, 



and so the first term on the right hand side of eqn. ( |3.35 ) can be dropped for small strain and 
eqn. (3.35) reduces to 



(3.39) 



/i3(B3 



^2^1 



10 



If Aj {i = 1,2,3,^ or no subscript) is the stretch, in one dimension, corresponding to the defor- 
mation gradient Fj, then Xf and Aj/Aj are the equivalent values in one dimension, corresponding 
to Bj and Dj. If e^ is the true strain for the stretch Aj, then e^ = InXi and so, q = Xi/Xi. Hence, 
eqn. (3.39) reduces to 



/i3 (A^ - 1) 



A, 



V2 



X, 



(3.40) 
or 

(3.41) /i3 (e'^^ - 1) = r/2e2, 

which under the assumption of small strain (i.e., es ^ 1), reduces to 

(3.42) 2/^363 = ^262. 



Following a similar analysis, in one dimension, eqn. (3.29);, becomes 
(3.43) 2/i3e3 + 2fipep = r]iei, 



and eqn. (3.29 )a reduces to 

(3.44) a = 2/1363 + 2fipep, 

where a is the one dimensional stress. In addition, eqn. (3.3), reduces to 

(3.45) A = A1A2A3, 
and so 

(3.46) e = ei + 62 + €3. 



Similarly, eqns. (3.4), (3.5) reduce to 
(3.47) e = 



Cp + ei, 



62 + 63. 



The eqns. (3.42-3.44), (3.46), (3.47) are, in fact, the equations obtained if we have the spring- 
dashpot arrangement shown in Fig. [3ta). 



We shall now show that these equations (i.e., eqns. (3.42-3.44), (3.46), (3.47)) reduce to the 
form of eqn. (1.1). Now, differentiating eqn. (3.47)b with respect to time and using eqn. (3.42), 
we obtain 



(3.48) 



. 2/i3 ^. 
ep = £3 + €3. 

V2 



Also, differentiating eqn. (3.44) with respect to time and dividing by 772, we find 
(3.49) - = -^€3 + 



2Ai3 . 2/i„ . 



V2 



V2 



V2 



~P1 



and differentiating eqn. (3.44) twice with respect to time and dividing by 2/^3 leads to 

a 



(3.50) 



2/i3 



£3 + 



fJ'3 

We add eqns. (3.49), (3.50) and use eqn. (3.48), to get 
(3.51) - + ^^ = Z£2ie; + ( 1 



-p- 



V2 2/X3 



V2 



jlp 
At3 



up. 
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From eqns. ( [3l3| ), ^A^ 

(3.52) 



a = rjiei = r]i[e- ep) , or, ep = e 



Vi 



Using eqn. (3.52) in eqn. (|3.51|) leads to 
2/ip 



(3.53) 



-a 



1 , /ip 



a 2fip . 



+ — ]^+7^ = ^^e + I 1 + ^ I e, 
^71^2 V^i l^sVi V2J 2/X3 r]2 V /^3. 



f^p\ . 



which can be re-written as 

^72 , ^72 , Vi 



(3.54) 



CT + I 1 1 a H a = rjie -\ 1 H — - ] e. 

2/ip 2/^3 2/ipy 4/ip/i3 2/ip V /U3, 



Eqn. (3.54) is in the same form as eqn. (1.1), with 



(3.55) pi = -— + -— + -— , p2 = -. , gi = Vu Q2 = ^ — 1 + — 

Zfip 2/i3 2fip 4/ip/i3 2fip \ fi3 



ni 



AAA 



11 



TI2 



2^3 



(a) 



12 -A/V^- 



(c) 



2^-3 



1G 



AAA^ 

2|l2 



1]3_ 2)Jl4 



Til 



■A/V^ 

2)i2 



(b) (d) 

Figure 3. Various spring-dashpot arrangements which reduce to the one- 



dimensional Burgers' fluid model (eqn. (1.1)). 



4. Model 2 

4.1. Preliminaries. Once again, let kh denote the reference configuration of the body. We 
shall assume that the body has two evolving natural configurations (denoted by t^piit), f^p2it)), 
but the manner in which they store the energy is different from that considered previously, with 
Fi,i = 1,2,3, being the gradients of the motion as discussed in model 1. We shall also use the 
definitions in eqns. (3.1), (3.2). Thus, eqn. ( 3.3[ ) applies here too. In addition, let us call the 
gradient of the motion from kr to Kp2{t) by Fc (see Fig. 111). It immediately follows that 

(4.1) Fg = F2F1. 



12 




>i(i) 



'P2(t) 



Figure 4. Schematic to illustrate the natural configurations for model 2. The 
body dissipates like a viscous fiuid during its motion from, kr to Hp^it), and kr 
to Kp^((). The body stores energy like a neo-Hookean solid during its motion from 

Kp^(t) to Kp2(t) and Kp^(t) to Kt. 



We shall denote the velocity gradient and its symmetric part corresponding to Fq by 



(4.2) 



Lg :— FgF^ , Dg :— - [Lg + ^g) ■ 



Also, from eqns. (3.3), (4.1), 

(4.3) F = FsFg. 

Following a procedure similar to the one followed previously for model 1, it can be shown that 



(4.4) Dg = D2 + -(FaLiF^-^ + F^'^Lf F 



Tt Tt:^T\ 



D = Ds + -(FsLgF:;! + F^^L^fT' 



'3 -r -l^raiJG-Ts 



-'G-^ 3^5 



along with 

(4.5) 



-2F,DGFr 



B 



'si-fGrs. 



B. 



-2F2D1F: 



2 5 



where A:= A — AL^ — LgA is the Oldroyd derivative when the natural configuration Kp^(^t) is 
the current configuration. In addition, from eqns. (4.4) and (4.5), along with the assumption 



that F2 = V2, F3 = V3 in virtue of the body being isotropic, we get 

(4.6) I.B2 = 2B2.(Dg-Di), I.B3 = 2B3.(D-Dg). 

These relations should suffice for our calculations for studying the response of model 2. 

4.2. Constitutive assumptions. In this model, we shall assume ip, and ^, to be of the form 

(4.7) 7/> = ^(B2,B3), ^ = e(Di,DG). 
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Now, assuming that the instantaneous elastic responses are isotropic and the body is incom- 
pressible, we choose 



(4.8) 

and 
(4.9) 



V^(B2, B3) = g(I.B2 - 3) + g(I.B3 - 3), 



e(Di, Dg) = r^iDi.Di + vg^g-'Dg. 



The above assumption implies that the body possesses instantaneous elastic response from the 
current configuration Kj to the natural configuration K,p2(t) and from the natural configuration 
Kpj(j) to the other natural configuration Kp^(ty It stores energy like a neo-Hookean solid during 
these two motions. In addition, the responses from the two natural configurations (Kpi(t), i^pi(t)) 
to the reference configuration Kn are purely dissipative, similar to a linear viscous fiuid. In fact, 
the response of the body as it moves from kr to K,p2{t) is similar to that of a "variant" of an 
Oldroyd-B fiuid (see Rajagopal and Srinivasa j^) i.e., the natural configuration /€p2(t) evolves 
like that of an Oldroyd-B fiuid with respect to the reference configuration kr. 



On substituting eqns. (4.8), (4.9) in (2.12), using eqn. (4.6) and simplifying, we get 
(4.10) (T - ^3B3).D + (/i3B3 - /i2B2).DG + /i2B2.Di = r/iDi.Di + vg^g-^g- 



Since, the right hand side of eqn. (4.10) does not depend on D, the incompressibility constraint, 
trCD) = 0, leads to 



(4.11) 



-pi + /igBa 



where —pi is the reaction stress due to the incompressibility constraint. Using, eqn. (4.11) in 



(4.10), we must have 



(4.12) 



{fi3^3 - yU2B2).DG + yU2B2.Di = T^iDi.Di + TjgBg-'Dg- 



Now, we maximize the rate of dissipation by varying T)i,T)g for fixed B2,B3 with the con- 
straints 

(4.13) tr(Di) = 0, tr(DG) = 0. 
Finally, we arrive at the following set of equations: 

T = -pI + /X3B3, 

(4.14) /X3B3 - /i2B2 = -p'l + Vg'Dg, 

/i2B2 = -/I + r/iDi, 

where p, p', p" are the Lagrange multipliers with 
(4.15) 



[/i3tr(B3) - /i2tr(B2)] 



V 



--/i2tr(B2) 
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Pre-multiplying and post-multiplying eqns. (4.14);,, (4.14 )c by V3 and V2 respectively, 



eqn. (4.14) reduces to 



(4.16) 



f^sBl - /iaVsBsVa = -p'B^ - ^ B3, 
fi^Bl = -p"B2 - I B2, 



with eqn. (4.15). If we denote H2B2, jJ^sB^ by 81,82 respectively, then the final constitutive 
relations for this model are 

T = -J9l + 82, 

(4.17) 



n 1 r)i ^G 

8? = -tr(808i-|8i, 



1 
3 



8^ - VS281 V82 = =: [ir(82) - tr(8i)] 82 - ^ 82 • 



2 



4.3. Reduction of the model to one dimensional Burgers' model. For simplicity, we shall 



use eqn. (4.14) for the reduction. Now, eqns. (4.14 )6^c can be re-written as 



(4.18) 



/i3 (B3 - I) - /i2 (B2 - I) = 3 [/^3 {tr{Bs) - 3) - /i2 (tr(B2) - 3)] I + r^cDc, 
//2 (B2 - I) = ^fi2 (tr(B2) - 3) I + r^iDi. 



Assuming that the displacement gradient associated with elastic response is small, leads to 

(4.19) ||B,-I|| = 0(7), 7<1, t = 2,3. 



The first term on the right hand sides of eqn. (4.18)fe^c can be neglected. Then, eqn. (4.18) reduces 
to 



(4.20) 



/i3 (B3 - I) - /i2 (B2 - I) = VG^tG, 

/i2 (B2 - I) = r^iDi. 



In one dimension, eqn. (4.20) becomes 
(4.21) 



A 



G 



/i3(A^-l)-/i2(A^-l)=r/Gv^, /i2(A^-l) 

■^G 



Vi 



Ai 
Ai 



where Aj {i = 1, 2, 3, G or no subscript) is the stretch, in one dimension, corresponding to the 
right stretch tensor Vj. Using InXi = e^ (e is the true strain), with the assumption of ej -C 1, 



eqn. (4.21) reduces to 



(4.22) 2;U3e3 - 2^^262 = VGe'c, 2/i2e2 = r/iCi. 
In addition, eqn. (4.14)^ reduces to 

(4.23) a = 2/i3e3. 



Eqn. (4.3) together with eqn. (4.1), in one dimension, reduces to 
(4.24) e = eG + e3, or, eG = e2 + ei. 
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dissipative response 

\Fi 



>i(t) 



Mt) 



Figure 5. Schematic to illustrate the natural configurations for model 3. The 
body's response is viscous fiuid-like and elastic solid-like, during its motion from, 
kr to Kpt(t), and i^p2{t) to kj respectively. From Kp^{t) to Kp^it), the response is 
Kelvin- Voigt solid-like. 



The spring-dashpot arrangement in Fig. 3[b) also yields eqns. (4.22), (4.23) along with 



eqns. (|3.46|), (|4.24|). These equations on simplification reduce to 



(4.25) 



a + 



2/i2 2/^3 2fl3 



a 



cr={r]i+ r]G) e + ^^ e, 



4/i2/i3 



2/i2 



which is same as the Burgers' one dimensional model (eqn. (1.1)), with 

ViVg 



(4.26) 



2/i2 2/i3 2/i3 



P2 



4/U2/i3 ' 



qi=Vl+ VG, 92 



VlVG 

2/i2 ■ 



5. Model 3 

5.1. Preliminaries. As with models 1 and 2, for this model we shall assume that the body has 
two evolving natural configurations {Kp-^(^t),i^pi{t), see Fig. [s]). We shall also use the definition of 



1,2,3 used for the previous models in addition to the definitions in eqns. (3.1), (3.2) and 



the relation eqn. (3.3). Further, we shall also choose F2 = V2 and F3 = V3. We recall from the 



preliminary discussion concerning models 1 and 2, that 



(5.1) 



D2 



-Vz'Bsy^', Di 



1 1 Vg -, 

-V2-I B2V2-\ 



These definitions and relations shall be used in the following analysis. 

5.2. Constitutive assumptions. For this model, we shall assume the specific stored energy, ip 
and the rate of dissipation, ^ to be of the form 



(5.2) 



^^V^(B2,B3), e = e(Di,D2 
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Specifically, in virtue of the body being incompressible and isotropic, we choose, 

(5.3) ^(B2, Bs) = ^(I.B2 - 3) + f (I.B3 - 3), 

ip Ip 

and 

(5.4) e(Di, D2) = ViDi.Di + r/p2.D2, 

i.e., the body possesses instantaneous elastic response from the current configuration Kt to the 
natural configuration i^p^{t) and stores energy like a neo-Hookean solid. Also, the response of the 
body between Kpi(i) to ^^■p2{t) is similar to that of a Kelvin- Voigt solid. The body also dissipates 
like a linear viscous fluid during its motion from k^ to i^p^(t)- 



On substituting eqn. (5.3) into eqn. (2.12) and using eqn. (3.18) we get 



(5.5) T.D - /isBa.Ds - /isBs. 
which reduces to 



D - D2 - - (FaLiF^-^ + F^-^Lf F^) 



/"3, 



T^iDi.Di+r/^Ds.Da, 



(5.6) (T-;f/3B3).D + (/i3B3-/i2B2).D2 + ^B3. (F2LiF2-^ + F^^Lf F^) = rl^Vi^Ji^^rl^Ti^X^^. 



Using eqn. (5.6), we maximize the rate of dissipation with incompressibility as a constraint, i.e., 
(5.7) tr(D) = tr(Di) = trij^i) = 0, 

by varying D, Di, D2 for fixed B2, B3 and get: 

T = -pI + /i3B3, 
P3B3 - /i2B2 = -p'l + ?72D2, 

^ (F^BgFa"^ + F2-^B3F2) = -/I + r/iDi, 



(5. 



where p, p', p" are the Lagrange multipliers with 

, 1 



(5.9) 

and 
(5.10) 



-p 



~P 



[/i3tr(B3) - P2tr{B2)] 
/i3tr(B3), 



v'Ai 



/i3B3.F2WiF2-^ 

r/;Di.Di+ 77^2.02 



1,2. 



eqn. (5.8 )c by V2 and using eqn. (5.1), we find that 



Pre-multiplying and post-multiplying eqn. (5.8)?, by V3, pre-multiplying and post-multiplying 



'pl + ;U3B3 



(5.11) 



/i3B2 - /i2V3B2V3 = 


-p'B, - 


-yB3, 


^(B2B3 + B3B2) = 


-p"B, - 





along with eqn. (5.9). 
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If we call fi^B^, (12^2 by 81,82 respectively, then, the final form for the constitutive relation 
can be given as 

T = -j9l + 8i, 

1 



(5.12) 



1 



8182 VSi = 3 [tr(8i) - tr(82)] 81 - | s[, 



(8281 + 8182 



1 r]i ^G 

-tr(8i)82 - y 82 . 



5.3. Reduction of the model to one dimensional Burgers' model. Following the method 



used in 4.3, eqn. (5.8), in one dimension, reduces to 

a = 2/1363, 

(5.13) 2/1363 - 2/i2e2 = 77262, 

2/i3e3 = ^iei- 

The above set of equations, can also be obtained from the spring-dashpot arrangement in 
Fig.gc). 



Now, eqn. (5.13) can be re-written as 
(5.14) 



a = 2/i3e3, a = 2/i2e2 + '72e2, cr = r^iei. 



Also, differentiating eqn. (3.46) with respect to time and using eqns. (5.14)a^c, we obtain 

0" 



(5.15) 



a 
Vi 2/i3 



Now, multiplying eqn. (5.15) with 2/i2, multiplying the derivative of eqn. (5.15) with respect to 



time with r/2; then, adding these two equations, along with eqn. (5.14)?,, we get 

(5.16) a +H \ ] cr + - — a = 2/i2e + r]2e, 

Vi \ ^71/^3/ 2/13 

re-written as 
(5.17) 



_2/i2 2/i2 2/13/ 4/i2//3 



Vim .. . , ??1??2.. 

a = r]ie+ - — e. 



2/12 



Thus, eqn. (5.17) has the same form as eqn. (1.1), with 

Vi 



(5.1^ 



2/12 2/12 2/13 



P2 



^1^2 



4/i2/i3 ' 

6. Model 4 



qi = vu 92 



^1^2 
2/12' 



6.1. Preliminaries. Once again, we shall assume that the body has two natural configurations 
associated with it, denoted by i^p^{t),f^p2{t)- However, in this model, the evolution equations of 
the two natural configurations are not coupled and they evolve independently (see Fig. |6]). We 
shall denote the gradients of the motion from kr to Kpi(t) and from tip^[t) to kj by Fi,F2. We 
shall also denote the gradients of the motion from kr to Kp^[t) and from Kp^(t) to Kj by F3, F4. It 
follows that 



(6.1) 



F2F1 = F4F 



4-^3- 
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Figure 6. Schematic to illustrate the natural configurations for model 4. The 
body's response is similar to that of a "mixture" of two Maxwell-like fiuids with 
different relaxation times. 



The left stretch tensor, velocity gradient and its corresponding symmetric part are denoted by 

1 



(6.2) 



B. := F.Ff , 



L, := F.F 



-1 



D, : = 



Also, a straightforward calculation leads to 



(6.3) 



V 
B2 = 



-2F2D1F: 
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V 
B4 = 



(L. + Lf) 



-2F4D3F[. 



1,2,3,4. 



^(B2, B4) = g(I.B2 - 3) + g(I.B4 - 3), 



6.2. Constitutive assumptions. Here, we shall assume the specific stored energy, ip and the 
rate of dissipation, ^ to be of the form 

(6.4) ij = ij{B,,B,), e = e(Di,D3). 

As the material is isotropic and incompressible, we choose, 

(6.5) 

and 

(6.6) e(Di, D3) = ^iDi.Di + r^aDs.Ds. 

This means that the response of the natural configurations {Kp-^(^t) ■, >^p2{t)) from the current con- 
figuration is like that of a neo-Hookean solid and the response from the reference configuration 
to the natural configurations is similar to that of a linear viscous fiuid. Thus, Burgers' fiuid can 
also be perceived as a "mixture" of two Maxwell-like fiuids with different relaxation times. 
We shall set 

(6.7) F2 = V2, F4 = V4, 
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where V2, V4 are the right stretch tensors in the polar decomposition of F2,F4, based on the 



assumption of isotropic elastic response. Hence, from eqn. (6.3) and eqn. (6.7), we have 
(6.8) I.B2 = 2B2.(D-Di), 1.64 = 264.(0-03). 



On entering eqns. (6.5), (6.6) in eqn. (2.12) and using eqn. (6.8), we get 
(6.9) 



(T - /i2B2 - /i4B4).D + /i2B2.Di + /i4B4.D3 

Using the constraint of incompressibility 

(6.10) tr(D) = tr(Di) = triBs) = 



r/iDi.Di +77303.03. 



and eqn. (6. 9), we maximize the rate of dissipation by varying D, Di, D3 for fixed B2, B4 and get: 

T=-pI + I22B2 + /i4B4, 

(6.11) 1J2B2 = -p'l + r]iT>u 

Ai4B4 = -p"l + r]3^3, 
where p, p', p" are the Lagrange multipliers with 

, 1 



(6.12) 



p 



-/i2tr(B2 



1 



-p 



fj,4tr{B4^). 



Pre-multiplying and post-multiplying, eqn. (6.11);, by V2, and eqn. (6.11 )c by V4; then, using 



eqn. (6.3), we arrive at 



(6.13) 



/^2B^ 



l^4,B^ 



-pB2 - 

-/B4 



YB2, 
YB4, 



Eqns. (6.13 )a,6 represents the evolution equations of the natural configurations {i^pi{t),i^p2{t) re- 
spectively) . If we denote /i2B2,/X4B4 by Si, S2 respectively, then the final constitutive relations, 
for model 4, are 

T = -j9l + Si + S2, 



(6.14) 



S? = -tr(Si)Si 



S^ 



ySi, 



1 r7i V 

-tr(S2)S2 - I S2 



This model is a variation of the model proposed by Murali Krishnan and Rajagopal |8j. They 
considered stretch dependent dissipation, in constrast to our linear viscous fluid type dissipation. 

6.3. Reduction of the model to the one dimensional Burgers' model. For this model, 
we shall once again assume that the displacement gradient associated with the elastic response 
is small, and thus 



(6.15) 



Then, eqn. (6.11) becomes 
(6.16) 



IB, -II 



fJ'2(B2 
/i4(B4 



0(7), 7<1, 



2,4. 



T= -pl + /i2B2 + /i4B4, 

I) = mDi, 
I) = r/3D3, 
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which in one dimension reduces to 



(6.17) 



a = 2/1262 + 2/1464, 
2/1262 = r/iei, 
2/1464 = r/sea. 



Further, eqn (6.1), in one dimension, reduces to 
(6.18) ( 



62 + ei = 63 + 64. 



In fact, the spring-dashpot arrangement Fig. k3[d) leads to eqn. (6.17), (6.18). We shall now show 



that these two equations, on simplification lead to eqn. (1.1). Differentiating eqn. (6.18) with 



respect to time and using eqn. (6.17)b,c, we have 
(6.19) 



. 2/12 
6 = 62 + 62, 

Vi 

. 2/l4 . 

6 = 64 + 64. 

^3 



Eliminating 64 from eqn. (6.17|)a and eqn. (|6.19|)fe leads to 
(6.20) 



O' , O- 2/i2 /12 . 

6 = h 62 62. 

V3 2/i4 r/3 /14 



Solving eqn. (6.19)a and eqn. (6.20) simultaneously, we get 



(6.21) 



^2 



£2 



1 + ^ 


6 — 






2/^2 / 


M2 _ 


-'A 




m \ 


^M4 


VsJ 




/ \ 

1 + 21 


M- 


a 


2/^4 



rn _ M2 

m M4 



Now, differentiating eqn. (6.21 )a with respect to time and equating it to eqn. (|6.21|){,, we get 

6.22 ^ + U~ + ^^~ K + ~A ^ = ^1 + ^3 e + - — 1 + — I 6. 

V2/i2 2/14/ 4/13/14 2/l2 V /^4. 



This is of the same form as eqn. (1.1) with 



(6.23) 



2/12 2/14 4/13/14 



gi = ^1 + ^3, g2 = 7^ — 1 + — 

2/12 V At4 



7. Final Remarks 



We have shown four sets of energy storage and rate of dissipation which lead to four different 
three dimensional constitutive relations, which reduce in one dimension to the model developed 



by Burgers (eqn. (1.1)). Each of these three dimensional models can claim equal status as 
representing the three dimensional generalization of Burgers' model. We have chosen two natural 
configurations instead of one in all of these models. This is to incorporate two relaxation times 
possessed by Burgers-like fiuid bodies. For example, in an asphalt concrete mixture (which has 
been shown to exhibit Burgers-like fiuid behaviour), the aggregate matrix has a small relaxation 
time whereas the asphalt mortar matrix has relatively larger relaxation time (see [S]) and the 
choice of two natural configurations seems natural. It is possible that several other choices for 
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the stored energy and the rate of dissipation could lead to the same one dimensional model due 
to Burgers. Interestingly, the structure of the three dimensional models that we have developed 
are quite distinct. 

8. Appendix 
Now, from eqn. (3.21) and eqn. (3.32) in [T7] 



(A.l) m^^ll^J^f 



aDi 2 I 9Li V9Li 



rp\ 



Hence, using eqn. (A.l) 



d d 

-B3. (F2LiF2"^ + F2" L^ F2 ) = — — -2B3.F2L1F2 ^ 



9Di ^ ^^ ^^ ' ' " 9Di 



(A.2) 



d 



-2Fn B.sFn .Li 



[1 

[2 

[3: 
[4; 

[5 
[6 

[7; 

[8 

[9 
[10 

[11 
[12 
[13 

[14 

[15 
[16 



(9Di 

F^BaF^-^ + F^-^Fs. 
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